Nonsingular parametric oscillators Darboux-related to the classical harmonic oscillator 
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Interesting nonsingular parametric oscillators which are Darboux related to the classical harmonic 
oscillator and have periodic dissipative/gain features are identified through a modified factorization 
method. The same method is applied to the upside-down (hyperbolic) 'oscillator' for which the 
obtained Darboux partners show transient underdamped features. 
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It is a common fact to factorize the classical har- 
monic oscillator equation u" + lu 2 u = through 

(li + -^J (li ~ Ru i) u = °> wnere Rui = u'x/ui is the 
log-derivative of one of the solutions. We will use the 
harmonic mode u\ = cosuJot in the following. On the 
other hand, the equation [l[ 



Wo(2tan 2 uj t + l)v = 



(1) 



the result of the reversed factorization 
(-si — Rui) (4; + Rui) V = 0. As well known, equa- 
tion ([T]) is the Darboux-transformed equation of the 
oscillator equation. Since the frequency parameter 
is time dependent, it appears that the v oscillator is 
parametric with a periodic chirping. Unfortunately, the 
linear independent modes of this oscillator 
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are periodically singular which casts doubts on their pos- 
sible applications. 

In this Letter, we show that more general factoriza- 
tions as discussed in which we call a/3 factorizations 
herein, lead to a new class of Darboux type partner oscil- 
lators of the harmonic oscillator, which while being para- 
metric are not singular. Since these factorizations can be 
applied to any second-order linear differential equation, 
we present them next using the common variable x and 
then return to the time variables when we will discuss 
the simplest possible applications related to the harmonic 
and hyperbolic oscillators. 
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By a/3 factorizations we mean factorizations of second 
order linear differential operators performed with the fol- 
lowing type of first-order operators: 



Br =a-\x)±±[3(x), 
B+ =a(x)£+p(x). 



(3) 



The minus sign in the factorization operator B~ produces 
only some sign changes in the rest of the mathematical 
formulas which are not important for the final results 
and therefore we will consider only the plus sign in the 
following. 



Then, we have: 



B + B~ = D z 



h ap H — 

a a 



D + {l3 2 +a(3 x ) , (4) 



where D = -4- and the subindex x also means the deriva- 

ax 

tive with respect to x. 



The reversed factorization reads: 



B~B + = D 2 
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Thus, let us start with an equation of the form 

D 2 y + f( x )Dy + , g(x)y = , (6) 

which for definiteness we assume to have only regular sin- 
gular points in its coefficients. One can think of it to be a 
factorized product of either the first or the second form. 
However, in previous works 0, H| it has been found that 
the second form is really tractable, while the first fac- 
torization leads to more complicated nonlinear equations 
for the factoring coefficients. Thus, considering © to be 
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of the B B + type, one is led to the following system of 
coupled equations: 



a x + a 2 /3 + P = af, 

(3 2 + ^ = g- 



(7) 
(8) 



This system of equations can be decoupled and leads to 
the following Riccati equation 



z x - f(x)z + z 2 + g(x) = 



(9) 



in the dependent variable z — (5/ a. Thus, with one so- 
lution h(x) of this Riccati equation at hand, one can get 
ft = h(x)a and then ([7]) turns into the following cubic 
Bernoulli equation 



a x + h(x)a 3 + [h(x) - f(x)]a = 



(10) 



The solution of this equation is related to the solution of 
the first-order differential equation 



w x - 2(h ~ f)w -2h = 
through a = w -1 / 2 , which leads to 

e -r(h-f)dx' 



atf,h( x ) 
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(11) 
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where A is an integration constant. Then /3f t h(%) 
h(x)atf t h(x) and the reversed-factorized equation is 



D 2 y- 



~- + f(x) ) Dy+(g(x) + /3 x (a- a' 1 )) y = 0. 



(13) 

This equation is the Darboux-transformed partner of 
equation (j6|) and it will be of special interest in the fol- 
lowing since it will be used to introduce a one-parameter 
family of nonsingular parametric oscillators related to the 
trigonometric and hyperbolic oscillators. 

We move now to the simple but common case f(x) = 0, 
which means equations originally without the first deriva- 
tive or from which the first derivative has been eliminated 
according to a well known procedure. This is the stan- 
dard case in quantum mechanics and has been addressed 
in the particular case of the quantum harmonic oscilla- 
tor in @, Q . There is a slight difference produced by the 
minus sign taken in front of a in the operator B + which 
is included to get the negative sign in front of the ki- 
netic term in the one-dimensional Schroedinger equation 
and also a scaling factor of 1/ y/2 of the operators to get 
the 1/2 numerical factor in front of the second derivative 
which is present for physical reasons. The overall effect 
is a change of sign of the linear term in the Bernoulli 
equation (fTQ|): 



+ h(x)a 3 — h(x)a = , 



(14) 



with the solution easily obtained from (TT2"j) . The quan- 
tum harmonic oscillator case corresponds to h(x) = x Q 



if the particular solution is used, but one can also employ 
the general Riccati solution Q. 

We focus now on the harmonic oscillator equation in 
classical mechanics, for which in addition to f(t) = 
one has g(t) = const, and we change the independent 
variable from i to 1 We have two non-trivial cases: 

1. g(t) = u! 2 . This choice leads to the normal 
classical harmonic oscillator. In this case, we have [l[ 



hit) = wq tan uj^t , 



(15) 



where the initial condition of the Riccati equation is as- 
sumed to be included in the phase of the tangent func- 
tion. Then 



a = ± 



cos uj^t 



\/\ — cos' 2 coot 



ujq sinwoi 
\f\ — cos 2 uj t 



(16) 



We use the plus signs in the following since the minus 
ones do not produce any change in the results. 

The corresponding B + B~y = equation, i.e., y" — 



■13' (a 



L )) y — 0, takes the form 



y" + 2( o (t)ij y' + u; 2 (t)y = . 
In the coefficient of the first derivative, we identify 

Atanwo^ 
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A — cos 2 ujot 



(17) 
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as the so-called damping ratio [4|, which is a function of 
time, differently from the common damped oscillators for 
which it is just a number in the interval [0,1]. On the 
other hand, 



u 2 (t) 
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(19) 



is the time periodic frequency of this parametric oscilla- 
tor. The solution of (| 1 T[) is: 



y(t) = Ci 



coot + i sin 2wo£ 

2\/\ — COS 2 LjQt 
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2^/X — cos 2 toot 



(20) 



Notice that this solution does not have singularities if 
A ^ [0,1]. Another interesting fact is that if one chooses 
wq/2, one can write this solution in 



C 1 = _L an d C 2 
the form 



y(t) 
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Vcos 2 LJ{)t — A 



( Vl (t) + iv 2 (t)) 



(21) 



which shows the connections with the v oscillators. Plots 
of the real and imaginary parts of (|20|) for A = 2 are 
given in Fig. Q]and of its modulus in Fig. [21 while ( (t) 
and w 2 (£) are displayed in Fig. [3] and Fig. @] respec- 
tively. It is worth mentioning the fact that the friction 
coefficient in ()1T[) is negative in a periodic manner. In 
general, the friction coefficient can be negative only as a 
result of some energy pumping mechanism. Such cases 
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are known for example in plasma physics, especially for 
dusty plasmas @, but one can think also of negative 
differential resistance circuits as well as mesoscopic and 
nanoscopic transport processes. 

2. g{t) = — &q. This case corresponds to the upside- 
down harmonic oscillator in the following sense. The 
Riccati solution of ([§]) is 



h(t) = — fcotanh/co^ • 
For the standard factorization, one gets 

u" - /c 2 (2tanh 2 k t - l)u = 



(22) 



(23) 



as initial equation having the linearly independent solu- 
tions 
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coshfco^ ~ kocoshkot 

The Darboux partner equation is 
w — k^w = 
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(24) 



(25) 



with the well-known linear independent zero-modes w± 
cosh kot and W2 ~ sinh kot. 

The a/3 factorization coefficients are now 



a = ±- 



COsh knt 



ko sinh kot 



(26) 



(27) 



(28) 



cosh 2 k t - A 4(cosh 2 k t - A) 2 ' 

The general solution has the following explicit form 

= [(C 3 + iC 4 7r) + C 4 (2fc t + sinh2A: t)] 
4v A — cosh 2 kot 



\j A — cosh 2 fco^ \/A — cosh 2 fcoi 

The B + B~y = equation reads 

y" + 2a(t)fc y' + fc 2 (t)y = , 

where the damping ratio is now 

t . . Atanhfcnt 

Ch(t) = n — 

cosh kot — A 

and the parametric frequency-like coefficient is 
fc2 , 2 / I sinh2fc i 



If A < 1, no singularities occur in (|26)) - (p0| . Plots of the 
real and imaginary parts of the solution (|30p are shown 
in Figs. §5§ and (|6]), respectively, and of its modulus in 
Fig. 0. We also plot the damping ratio, (h(t), in Fig. flSJ 
and the frequency- like coefficient fc 2 (i) in Fig. (|9]), all of 
them for A = 0.5 and fco = 1. Examining these plots, 
it can be seen that after an initial period when both 
quantities are positive, the damping ratio drops to zero, 
while fc 2 (i) tends to — feg. In other words, the parametric 
equation (|27p turns into the standard hyperbolic equa- 
tion flu)) . 

We finally notice that if we write the B + B~ equation 

as 



D 2 y + F(x)Dy + G{x)y = 



then 



F = f-2- 



G = g + (3 X 



(31) 



(32) 



Thus, if a — ±1 as in the standard factorization no new 
equation is obtained. On the other hand, it is easy to 
check that for / = the case a = ±i leads to the un- 
damped but singular parametric oscillators obtained by 
the standard factorization for g = ± const., which also 
correspond to taking A = 0. More generally, in the frame- 
work of the a/3 factorization, one can get undamped para- 
metric equations only if / = 2—, which leads to F = 



and G 



1 



In summary, we have presented here the general 
method of a/3 factorizations and based on it we intro- 
duced an interesting family of nonsingular parametric 
oscillators which are Darboux related to the classical har- 
monic oscillators. The scheme has been also applied to 
the hyperbolic case. 
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FIG. 2: (Color online) Three-dimensional plot of the modulus 
of the solution ()20|) as a function of its real and imaginary 
parts for the same values of the parameters. 
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FIG. 3: (Color online) The damping ratio f (i) in (fT8l) for 
A = 2 and ujo = 3.5. 
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FIG. 4: (Color online) Plot of wg(t) given by equation (|19[) 
for A = 2 and ojo = 3.5. 
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FIG. 5: (Color online) The real part of the solution (|30[) for 
the superposition constants Cs — 2 and Ca = — 1 and ko — 1 
and A = 0.5. 
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FIG. 6: (Color online) The imaginary part of the solution 
(|30l l for the same values of the parameters. 
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FIG. 7: (Color online) Three-dimensional plot of the modulus 
of the solution (|30|) as a function of its real and imaginary 
parts for the same values of the parameters. 
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ko = 1 and A = 0.5. 




FIG. 9: (Color online) Plot of kg(t) given by equation 
for A = 0.5 and ko = 1. 



